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Abstract 

Research on combustion Is being conducted at 
Lewis Research Center to provide Improved analyti- 
cal models of the complex flew and chemical reac- 
tion processes which occur In the combustor of gas 
turbine engines and other aeropropulslon systems. 
The objective of the research Is to obtain a bet- 
ter understanding of the various physical proces- 
ses that occur In the gas turbine combustor In 
order to develop models and numerical codes which 
can accurately describe these processes. Activi- 
ties Include In-house research projects, univer- 
sity grants, and Industry contracts and are 
classified under the subject areas of advanced 
numerics, fuel sprays, fluid mixing, and radiation- 
chemistry. Results are highlighted from several 
projects. 

Introduction 

As we enter Into the age of the computer and 
the laser, new opportunities are emerging to make 
significant advances In the understanding of the 
physics of combustion. These advances are needed 
in the gas turbine engine Industry. Development 
costs and maintenance costs are prevalent factors 
when considering the next generation of gas tur- 
bine engines. The costs associated with the 
design and development of engine components. 
Including the combustor, are alarmingly high using 
today's design methodology.^ Also, the ever- 
increasing operating pressures and temperatures of 
new engines for Increased performance and better 
fuel economy are having adverse effects on dura- 
bility of hot section components. This results In 
high maintenance costs. Clearly what Is needed is 
an improved design methodology which can allow the 
design engineer to converge to the optimum design 
in a much shorter development cycle. 

At NASA Lewis Research Center, combustion 
research is integrated into activity in Internal 
fluid mechanics. The objectives of the research 
are to advance the understanding of flow physics, 
heat transfer and combustion processes which are 
fundamental to aeropropulslon, and to translate 
this knowledge into models and numerical codes of 
aerothermodynamic phenomena. The overall goal is 
to bring internal coiTiputational fluid mechanics to 
a state of practical application to propulsion 
systems. These models and numerical codes would 
then be available to the Industry to Incorporate 
Into their own engine/component design systems. 

The approach at Lewis Is to establish an 
Integrated computational-experimental research 
program. The activity consists of research on 


numerical methods, well-defined experiments for 
code and model verification, and the demonstration 
of computational codes for propulsion system 
components. In the area of combustion, research 
has been underway for some time In these three 
areas. with a recent reorganization In the 
Aeronautics Program at Lewis, this research activ- 
ity has been consolidated with similar research In 
other components such as compressors, turbines, and 
transition ducts and Is continuing In the Internal 
fluid mechanics research program. Although the 
amount of work currently being supported In com- 
bustion has been decreased In the new organiza- 
tion, the prime research needs In combustion are 
recognized and the program Is focusing on them. 

Figure 1 Is a cut-away view of a ficticious 
combustor which illustrates the complex fluid 
mechanics and combustion features. The flows are 
highly three-dimensional with turbulence levels in 
many cases comparable in magnitude to the bulk 
velocity. Liquid fuels are Injected as a spray 
which then undergoes vaporization and mixing. 
Chemical reaction occurs which causes changes In 
density and fluid mechanics properties and can 
cause the formation of a solid phase (soot) with 
its attendant high radiation heat transfer. The 
understanding of these physical processes is 
needed before accurate numerical codes can be 
built and used as a predictive tool in the design 
process. In addition the numerical methods for 
three-dimensional flows need Improvements In 
accuracy and efficiency In order to properly simu- 
late the features of these flows. This then is 
the framework of the current combustion research 
prnoram. As shown in Fig. 2, the activities are 
focused in four areas: advanced numerics, fuel 

sprays, fluid mixing, and radiation-chemistry. 

Both experiments and computational research are 
being conducted, with the long range objective of 
providing the numerical codes that can be used 
with confidence as a predict'.* tool in the indus- 
try's design system. In this report a research 
project from each of these areas will now be high- 
lighted as examples of the kinds of work currently 
being supported at Lewis. 

Advanced Numerics: Error Reduct n Program 

A major restriction to the development of a 
computer based design methodology is the accuracy 
of the numerical methods used in combustor flow 
codes. The upwind differencing currently used in 
these flow codes, ^ Introduces an appreciable 
error in the calculated result. This error (or 
numerical diffusion) Is frequently of such a large 
magnitude that It obscures the turbulence model 
used In the calculation. 



To alleviate this problem, NASA has conducted 
a program to Identify and Incorporate an Improved 
accuracy differencing scheme Into a combustor flow 
code. Under a portion of this program a variety 
of differencing schemes were examined In several 
test calculations. The schemes examined Included 
QUICK (Quadratic Upstream Interpolation) and SUD 
(Skewed Upwind Dlfferendog). QUICK differencing 
was developed by Leonard. This scheme Improves 
the accuracy of convective differencing by per- 
forming an upwind biased quadratic Interpolation. 
This scheme Is second order accurate and can pro- 
duce nonphysical oscillations in the solution. 

SUD, developed by Ralthby,” attains high accu- 
racy by differencing In an upwind manner along the 
flow streamlines. While maintaining the same for- 
mal accuracy as upwind differencing the truncation 
error In SUD Is smaller In magnitude. As with 
quick, sue can produce nonphysical oscillations In 
the so1i;t1on, therefore, a scheme to "bound* SUD 
was also examined. This scheme employs tho con- 
cept of flux-blending,' wherein a bounded flux 
determined from upwind differencing Is blended 
with the unbounded - but more accurate - SUD flux. 
The main factor Is to blend as little of the lesser 
accurate scheme while still maintaining a properly 
"bounded" solution. This procedure called BSUDS, 
starts from an Initial, totally skew differenced 
estimate and blends an upwind flux If the solution 
Is out of the range of neighboring values. If the 
solution Is In range, I.e., bounded, then no blend- 
ing Is performed. 

An Illustration of the accuracy of the 
upwind, QUICK, and SUD schemes Is seen In Fig. 3. 
The figure displays the results of a single point 
scalar transport calculation made for various flow 
angles. All schemes agree with an exact solution 
(no error) at a flow angle of zero, but departing 
from this each scheme displays some degree of 
error relating to numerical diffusion. The error 
displayed by upwind dlfferenclng^lncreases with 
flow angle up to a maximum at 45*. QUICK displays 
a similar behavior but the overall level of error 
IS much less. SUD displays an error maximum 
around 15* but tends to zero at angles approaching 
45 . Both QUICK and SUD display a much higher 
level of accuracy than upwind. 

A scalar transport calculation Is useful for 
generally examining some aspects of differencing 
scheme performance, but a more complete test Is a 
laminar flow calculation. The results of a series 
of laminar flow calculations from Ref. 8 are dis- 
played In Fig. 4. In this figure axial velocity 
profiles at a distance of one-half a duct height 
from the Inlet are shown for two different compu- 
tational meshes. In these calculations an Indica- 
tor of accuracy Is the steepness of the velocity 
profile. Steep velocity profiles are exhibited by 
QUICK and BSUDS with the upwind profiles exhibit- 
ing a high degree of numerical diffusion. On the 
coarse mesh, BSUDS appears to be more accurate 
than QUICK, while :r. a fine mesh there Is not much 
of a distinction between the two schemes. 

An Important aspect of Improved accuracy Is 
the effect that the differencing scheme has on the 
rate of convergence. The computational times to 
converge the system of governing equations In the 
previous laminar calculations Is shown In Table 1. 
The convergence times are ratloed to the upwind 
convergence times to clearly Illustrate the compu- 
tational penalty paid to attain Improved accuracy. 


Generally, the Improved accuracy schemes required 
from 3 to 15 times longer to reach a converged 
solution. To a decree one can hope that this 
coaiputatlonal penalty can be offset with the use 
of coarse meshes to achieve the same overall level 
of accuracy. In reality It appears tnat a rela- 
tively fine mesh Is needed even with the high 
accuracy schemes. In any case the need for 
Improved solution algorithms with these more accu- 
rate differencing schemes Is strongly demonstrated. 

Fuel Sprays; Study of Dilute Turbulent 
Particle-laden Jets 

A research topic of current Interest In com- 
putational fluid mechanics Is spray processes In 
combustion chambers. However, due to the compli- 
cations Involved In the spray processes, e.g., 
hydrodynamics of spray atomization, polydl sparse 
drop diameter distributions, drop shattering and 
coalescence, and evaporation and combustion of 
drops, spray processes In combustion chambers are 
excessively complex for systematic model develop- 
ment and evaluation. In contrast. Injection and 
mixing of monodlsperse solid particles Is a much 
simpler process, and facilitates the modeling of 
the underlining physics Involved In practical mul- 
tiphase flows. 

The objective of this work was to evaluate 
models of dilute turbulent particle-laden jets. 
Dilute sol Id-particle-laden jets In a still envi- 
ronment were considered. These flows are rela- 
tively simple and only Involve Interphase momentum 
exchange. Modeling efforts concentrated on effects 
of turbulent fluctuations on momentum transfer 
between the phases, as well as the dispersion of 
the particles by turbulent fluctuations. Three 
typical models of the flow processes were consid- 
ered: (1) a locally homogeneous flow (LHF) model, 
where slip between the phases was neglected; (2) a 
deterministic separated flow (DSF) model, where 
slip was considered but effects of turbulence on 
particle motion were Ignored; and (3) a stochastic 
separated flow (SSF) mode! where effects of both 
Interphase slip and turbulence on particle motion 
were considered using random-sampling techniques. 

The LHF model represents the simplest treat- 
ment of a multiphase flow. In this model, the 
particles and the continuous phase are assumed to 
have equal rates of turbulent diffusion. This 
Implies that Interphase transport rates are Infi- 
nitely fast, so that both phases have the same 
velocity at each point In the flow. The LHF 
approximation Is valid only for flows containing 
small particles, where characteristic response 
times of particles are small In comparison to 
characteristic times of turbulent fluctuations. 
Faeth’ has Indicated that LHF models only ylelo 
accurate predictions for particle sizes smaller 
than most practical applications. Despite the 
limitations, LHF models have some important advan- 
tages. They require minimal Information concern- 
ing Initial conditions of particle properties, 
which are often difficult to obtain for practical 
flows. The theoretical model of the flow Is equiv- 
alent to that of a single-phase flow and effects of 
multiple phases only appear In the representation 
of state relationships, bypassing the difficulties 
Involved with modeling Interactions between the 
phases. 
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The DSF model considers slip between the 
phases but neglects turbulent particle dispersion. 
In this model particles follow deterministic tra- 
jectories and only Interact with mean properties 
of the continuous phase. The computation Involves 
dividing the particulate phase Into representative 
samples whose motion and transport are tracked 
through the flow field using a Lagranglan formula- 
tion. An Eulerlan formulation Is employed to 
solve the governing equations for the continuous 
phase. The effect of Interphase transport Is con- 
sidered by Inserting appropriate source terms In 
the governing equations for the continuous phase. 
The source terms are determined using the “Particle 
Source In Cell* technique. Neglecting turbulent 
dispersion Is appropriate for flows containing 
large particles, where particle response times are 
large In comparison to characteristic turbulent 
fluctuation times. 

Most practical particle-laden flows and 
sprays, however, exhibit properties between the 
two limits represented by LHF and DSF models and 
therefore require consideration of turbulent par- 
ticle dispersion. The stochastic separated-f low 
(SSF) approach, first proposed by Gosman and 
loannldes^^ and subsequently developed by Faeth 
and coworkers, was adopted to study the 
effects of turbulence on particle dispersion. The 
SSF model Involves finding trajectories of statis- 
tically significant samples of Individual parti- 
cles as they move away from the Injector and 
encounter a random distribution of turbulent 
eddies. Particle trajectory computations are 
similar to the DSF model, except that instanta- 
neous eddy properties repla'.e mean-gas properties. 
Eddy properties are obtained by constructing the 
probability density function (PDF) of velocity, 
assuming Isotropic turbulence, and random sampling 
to find each velocity component. 

The three models of particle-laden jets have 
been systematically evaluated and some typical 
results are shown here. Comparisons of model 
predictions with experimental data are shown in 
Mgs. 5 and 6, taken from Ref. 13. Figure 5 
presents radial profiles of mean and fluctuating 
properties for the gas phase, and Fig. 6 presents 
radial profiles for the particle phase. Since 
this jet has a relatively low particle loading, 
gas phase properties are relatively independent of 
model type and both the LHF and SSF models are In 
good agreement with the measurements. For the 
particle properties SSF predictions agree reason- 
ably well with the data. Predictions of the LHF 
model for par' icle properties are not at all 
satisfactory, since neglect of slip overestimates 
particle fluctuations levels and rates of spread. 
The DSF model also yields poor results because 
neglecting particle dispersion by turbulence 
causes the rate of spread of the particles to be 
substantially underestimated. Similar results 
were found throughout the data base used in the 
study. 

For the two separated-f low models used in 
this study (DSF and SSF) the model predictions are 
relatively insensitive to the specification of 
gas-phas9 initial conditions. The specification 
of initial particle properties, however, exerts 
much more pronounced effects on predictions. 
Accurate measurement of particle initial condi- 
tions, e.g., particle size, mean and fluctuating 
velocities and particle mass flux, is essential to 


obtain meaningful predictions of separated-f low 
models, especially when particle relaxation time 
Is large compared to particle residence time In 
the flow field. 

Practical particle-laden flows and sprays 
usually Involve regions where the dilute flow 
approximation Is not appropriate, such as near the 
Injector. Phenomena to be considered for future 
model extension In application Involving dense 
particulate flows Include: the finite volume 
fraction of the dispersed phase, effects of nearby 
particles on particle transport properties, and 
effects of particle collision and coalescence. 
There Is also a pressing need to overcome the 
experimental difficulties encountered In the dense 
flow region so that further progress can be made 
toward developing reliable m(v;*ls for this region. 

Fluid Mixing; Direct Numerical Simulations of 
Chemically-Reacting Mixing Layers 

The advent of high speed vector processing 
computers has made possible the study of turbu- 
lence and turbulence-chemistry interactions via 
computational "experiments." One promising 
approach that has been supported at NASA Lewls^^ 

Is Direct Numerical Simulation (DNS). Direct 
Numerical Simulations solve the complete time- 
dependent, Navier-Stokes equations at a range of 
flow Reynolds number wher.*1n all the scales of 
turbulence can be fully resolved on a computa- 
tional mesh. This typically requires the calcula- 
tion of flows where the Reynolds number based on 
the Taylor macroscale is of the order of 100. The 
key assumption involved is that the large scales 
of turbulence which can be represented on the com- 
putational mesh will not significantly vary with 
Reynolds number. While this may not be true for 
all the aspects of turbulent flow which one might 
desire to study, a previous study of Incompres- 
sible turbulent flows“ has indicated that a 
great deal of Insight Into the structure of turbu- 
lent flows can be obtained. One of the main goals 
Is that this information will be effective In 
developing turbulence-closure models. 

Incompressible turbulence models typically 
require Information on fluctuating pressure- 
velocity correlations that are difficult to 
measure experimentally. Reacting flows further 
complicate the problem with the need for an 
Instantaneous measure of density and the number of 
double and triple correlations rapidly becomes 
overwhelmino. Indeed, even the form of averaging 
(Fauvre or Reynolds) to be used in constructing 
the governing equations Is not straight-foward.^” 
The promise of DNS is that it can provide complete 
flow field values that can be processed to provide 
any necessary correlation or test various closure 
models. 

The above states the promise of DNS, but 
there are a number of ''•ealitles" which should 
also be mentioned. The first has already been 
noted: simulations must be made at low Reynolds 

number. This is primarily a numerical resolution 
limitation. Even employing highly accurate 
psuedo-spectral methods^' and large computa- 
tional meshes (typically 64 by 64 by 64 for the 
results reported herein) the range of turbulent 
scales that can be represented is limited. A 
second limitation is that the numerical methods 
used restrict the calculation to relatively simple 
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geometries. The numerical calculations discussed 
here are of a time-evolving shear laiyer as opposed 
to the experimental 1v measured special mixing 
layer. The time-evolving mixing layer Is not a 
direct analog of the special mixing layer, and 
some physical phenomenon cannot be represented. 

For example, the f1ov> streamlines In the time- 
evolving case are not aligned along the flow 
direction as they are In the special case. A 
final limitation relates to the statistical nature 
of DNS results. One simulation must typically run 
for several CPU hours to simulate a small fraction 
of physical time (on the order of a second). Sta- 
tistics sensitive to scale variation may not be 
aoequately sampled. It Is easily seen that more 
small scales of turbulence can be represented on a 
finite computational domain than larger scales and 
this can bias the statistical results. These fac- 
tors require care In the sampling of simulation 
results to obtain statistically stationary 
Information. 

UespUe the current limitations Inherent to 
DNS a number of turbulent flow phenomenon can be 
well represented. An example Is the experimental 
versus computational flow visualization results 
'hown In Fig. 7. The experimental flow visualiza- 
tion was provided by Koochesfahan1^° and dis- 
plays a flow structure resulting from the merging 
of two discrete vortices. This flow structure Is 
strikingly simlllar tu the results of a two- 
dimensional numerical simulation. The experiment 
was conducted at a Reynolds number and with fluids 
differing from the computation, yet the comparison 
Is encouraging. 

Another favorable aspect of DNS results Is 
displayed in Fig. 8. Similarity theory predicts 
that the mean velocity half-width of the shear 
layer should grow linearly. As seen In Fig. 8, 
the time evolution of the velocity half-width 
grows linearly with time In the numerical simula- 
tions. If a transformation relating space and 
time Is performed^' the time evolving results 
are found to be In good agreement with similarity 
theory. This agreement Is favorable for a large 
range of Initial conditions as shown on Fig. 8. 

A final comparison between simulated results 
ana experimental data Is shMn In Fig. 9. The 
experimental data, Mungal,^’ Is of a mildly exo- 
thermic chemical reaction. The similarity plot of 
Integrated product level from a numerical simula- 
tion agrees fairly well with the shape and peak of 
the experimental data. The shape of the product 
proflU Is somewhat sensitive to the Initial con- 
ditloni used In the simulation and current research 
Is being conducted to assess whether a1te'*nate 
initial conditions might clear up the small dis- 
crepancy displayed. 

The overall agreement between DNS results and 
experimental Information Is strongly encouraging. 
The highlights of the current research Is merely 
Illustrated here and Ref. 14 Is recommended for 
further details. 

From this army of favorable Initial results 
It appears that this approach can provide Insights 
Into the physics of turbulent chemistry which will 
one day lead to Improved closure models usable by 
combustor designers. 


Radiation, Chemistry; Fast Algorithms for 
bomputllon kinetics caicuiatiSnl 

The ordinary differential equations (DDE's) 
describing complex chemical reactions are charac- 
terized by widely different time constants. 

Although the DDE's are stable, standard numerical 
techniques such as the popular explicit Runge-Kutta 
and Adams methods are prohibitively expensive to 
use because of the severe steplength restrictions 
Imposed by the requirements for numerical stabll- 
1ty.‘® Such systems of DDE's are commonly 
referred to as "stiff systems. 

The aim of the present Investigation was to 
examine currently available Integration techniques 
for solving chemical kinetic rate equations. The 
motivation behind this work Is the Increasing 
Interest In (1) modeling the reaction mechanisms 
describing the consumption of fuels and pollutant 
formation and destruction, and (2) multidimensional 
modeling of reactive flows, which Includes the 
equations of fluid motion. The former results In 
the need to Integrate large systems of nonlinear 
DDE's. The latter results in the need to Inte- 
grate the DDE's at several thousand grid points. 

To make such calculations practicable. It Is 
necessary to have a very fast homogeneous batch 
chemistry Integrator. 

The techniques examined Include two general- 
purpose solvers, EPISDDE^^ ani' LSDDE,‘^ developed 
for an arbitrary system of DDE's, and the special- 
ized techniques CHEMEQ,^^ CREKlD,^® and GCKP84,2® 
developed specifically for chemical kinetics appli- 
cations. Testing of these codes was accomplished 
by application to two practical combustion kinetics 
problems. Both problems described adiabatic, homo- 
geneous, gas-phase, transient, batch combustion 
reactions at constant pressure. Test problem 1 
described the Ignition and subsequent combustion 
of a mixture of 33 percent carbon monoxide and 
67 percent hydrogen with IDO percent theoretical 
air. It consisted of 12 reactions among 11 spe- 
cies. Test problem 2, Involving 30 reactions among 
15 species, described the Ignition and subsequent 
combustion of a stoichiometric hydrogen-air mix- 
ture. Both problems were Integrated over a time 
period of 1 msec, which encompassed all three com- 
bustion regimes: Induction, heat release, and 
equilibration. 

In using LSODE, EPISODE, and CHEHEQ the tem- 
perature was calculated by employing two different 
methous. In method A the temperature was evalu- 
ated via an Iterative solution of the algebraic 
enthalpy conservation equation. In method 8 the 
temperature was computed by solving Its ODE. Both 
CREKID and GCKP84 were written explicitly for non- 
Isothermal chemical reactions and therefore Include 
calculation procedures for the temperature. 

Figures ID and 11 present the variation of 
the local error tolerance with the computer (I.e., 
CPU) time required to solve te*.t problems 1 and 2, 
respectively. In these figures CPU Is the com- 
puter time (In seconds) required on the NASA Lewis 
Research Center's IBM 37D/3D33 computer. All 
results were obtained using single-precision 
accuracy, except GCKP84 which was In double- 
precision. The error control performed Is dlff- r- 
ent for the different codes. Consequently, the 
local error tolerance, EPS, ooes not have the same 
meaning for all the codes. For LSDDE, GCKP84, 
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CREKIO. and CHENEO, EPS Is the local relative 
error tolerance/^ For EPISODE, however, EPS Is 
a mixed relative and absolute error tolerance - 
relative for species with Initially nonzero .X)le 
fractions (I.e., reactants) and for the tempera- 
ture, and absolute for species with Initially zero 
mole fractions (I.e., Intermediate species and 
products). 

Figures 10 and 11 show that LSODE Is the 
fastest code currently available for solving com- 
bustion kinetic rate equations. Because GCKP84 
has been designed for a variety of chemical kinet- 
ics problems It requires more work per step than 
LSOOE and Is hence slower. CREKlD and, especially 
for test prob’em 2, EPISODE are attractive alter- 
natives to LSODE. However, there are many prob- 
lems associated with using EPISODE for solving 
combustion kinetics rate equations - details are 
given In Ref. 20. Figures 10 and 11 show also 
that the use of temperature method A does not 
result in significant 1neff1c1e''cy. On the con- 
trary this method can be significantly faster than 
evaluating the temperature by Integrating Its ODE 
(Fig. 11). 

For multidimensional chemically reacting flow 
calculations computational speed Is of primary 
concern and moderate accuracy Is adequate. How- 
ever, In developing and validating reaction mecha- 
nisms, accuracy Is of c'ltlcal importance. An 
accuracy comparison of the techniques examined In 
this study was made as follows. For each test 
problem a mean Integrated root mean square error 
(Erms) MAS estimated. Ep^j Is approximately 
the average true error Incurred by a technique In 
solving the complete problem. Figures 12 and 13 
present the variation of Ep^j with the user- 
supplied value for the local error tolerance, 

EPS. Note the significant discrepancies between 
values specified for EPS and the errors actually 
obtained. For both probleins temperature method A 
Is as accurate as temperature method B - In many 
cases. It Is significantly more accurate. LSODE 
Is the most accurate code for solving chemical 
kinetic rate equations. However, 6CKP84, CREKlD, 
and ChEMEQ-A compare favorably with LSODE for 
small values of EPS. EPISODE is significantly 
less accurate than the other codes because the 
error control used In It Is Inappropriate for 
chemical kinetics applications.^'^ 

The present study has shown that LSODE Is the 
fastest and most accurate code currently available 
for solving chemical kinetic rate equations. How- 
ever, It Is recognized by combustion device mod- 
elers that LSODE Is not fast enough for economical 
calculations of multidimensional reacting flows. 
Current research efforts at NASA Lewis are aimed 
at further speeding up codes for solving chemical 
kinetic rate equations. 

Concluding Remarks 

The four research activities highlighted 
above are examples of NASA Lewis sponsored work 
which are building blocks toward achieving compu- 
tational codes of combustion and fluid mechanics 
which can be applied by Industry as predictive 
tools in their design system. Much work remains 
to be dene In this area and NASA Lewis will con- 
tinue to comriit itself to using Its resources In 
working with other segments of the technical com- 
munity In achieving these research objectives. 
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TABLE I. - RATIO OF CONVERGENCE TIMES 
FOR VARIOUS DIFFERENCING SCHEMES 
WITH UPWIND CONVERGENCE TIMES 
USED AS THE STANDARD 
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• HIGH TURBULENCE LEVELS * 2 PHASE WITH VAPORIZATION 

Figure 1. - Cutaway Illustration of a portion of a full annular 
combustor system, with major features highlighted. 
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Figure Z - Combustion research has long range objectives to produce 
predictive numerical codes that have practical application in engine 
manuticturer's combustion design system. 






Figure 3. - Results of a scalar transport test 
calculation for various differencing schemes. 
Error from an exact solution is displayed 
versus flow angle. 








Figure 5. • Radial variation of gas-phase 
mean and turbulent quantities in case 
3 particle-laden jet x/d - 40. SMD ■ 
119 urn. LR-a66. 
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Fk.ure 6. * Radial variation of solid-phase 
mein and turbulent quantities in case 3 
particle-laden jet x/d ■ 40. SMD - 119 pm, 
LR-a66. 
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(a) Comoutational results. tb) Experimental flow visualization. 


Figure 7. - Flow structure comparison between direct numerical 
simulation and experimental flow visualization. 
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Figure 8. - Mean velocity halfiwidth. (■ Zi/2>- 
versus time - three-dimensional simulations. 
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Figure 9. - Similarity plot of integrated product 
level versus radial location. 3-D simulation. 
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Figure IGL - Variation of the CPU time(s) with error 
tolerance. EPS, for test problem 1. All runs on 
IBM 370/3033. 
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Figure 11. - Variation of the CPU timeU) with error 
tolerance. EPS, for test problem 2. All runs on 
IBM 370/3033. 


AAEAN INTEGRATED ERROR 





loO 



LOCAL ERROR TOLERANCE, EPS 

Figure 13i. > Variation of the mean integrated global 
error (Ep,p Jwith the local error tolerance (EPS) 
for test problem Z 



